Combinatorics of Feynman Diagrams for the Problems with 

Gaussian Random Field 



E.Z.Kuchinskii, M.V.Sadovskii 

Institute for Electrophysics, 
Russian Academy of Sciences, Ural Branch, 
Ekaterinburg 62004-9, Russia 
E-mail: sadovski@ief.intec.ru 

Submitted to JETP, May 1997 



Abstract 



The algorithm to calculate the generating function for the number of 
"skeleton" diagrams for the irreducible self-energy and vertex parts is derived 
for the problems with Gaussian random fields. We find an exact recurrence 
relation determining the number of diagrams for any given order of perturba- 
tion theory, as well as its asymptotics for the large order limit. These results 
are applied to the analysis of the problem of an electron in the Gaussian ran- 
dom field with the "white-noise" correlation function. Assuming the equality 
of all "skeleton" diagrams for the self-energy part in the given order of pertur- 
bation theory, we construct the closed integral equation for the one-particle 
Green's function, with its kernel defined by the previously introduced gener- 
ating function. Our analysis demonstrate that this approximation gives the 
qualitatively correct form of the localized states "tail" in the density of states 
in the region of negative energies and is apparently quite satisfactory in the 
most interesting region of strong scattering close to the former band-edge, 
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where we can derive the asymptotics of the Green's function and density of 

states in the limit of very strong scattering. 

PACS numbers: 02.10.Eb, ll.10.Jj , 71.55.Jv, 72.15.Rn, 
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I. INTRODUCTION 

Different methods of summation of Feynman diagrams are widely used in rather wide 
class of problems of theoretical physics dealing with elementary excitations propagation in 
systems with static random fields, created by different types of inhomogeneites. A simplest 
example is the problem of an electron in a metal with impurities. Apparently this was the 
first problem in which the diagram technique considered in this paper was formulated for 
the first time. P]|| . Analogous formalism is being widely used in problems of statistical ra- 
diophysics and optics, dealing with the propagation of electromagnetic waves in disordered 
media Equivalent mathematical approach is used in a number of problems of critical 
phenomena in disordered systems [f§], in the problem of polymer chain with "excluded vol- 
ume" and in some other problems of polymer physics ||. The same diagram technique 
describes the regular model of critical phenomena with "zero-component" order parameter 

I- 

In any problem, dealing with summation of Feynman diagrams, any kind of information 
on combinatorics of graphs, i.e. on the number of diagrams of different types for the given 
order of perturbation theory, is quite useful. In this paper we present a detailed study of 
Feynman graphs combinatorics for the problems with the Gaussian random field. 



II. GENERATING FUNCTION FOR THE NUMBER OF "SKELETON" 
DIAGRAMS. RECURRENCE RELATION. 

For definiteness we shall always discuss the problem of an electron with energy E and 
momentum p, propagating in the Gaussian random field (random impurity system) flH0. 
The averaged Green's function is determined by the perturbation series shown in Fig. 1(a). 
In a standard way this series is reduced to Dyson's equation: 
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where e p = ^— - is the free electron spectrum, while the self-energy part T,(E,p) is defined 
by the "skeleton" graphs of Fig. 1(b), where the internal electronic line represents the fully 
renormalized ("dressed") Green's function G(E,p). 

The total number of diagrams in the iV-th order of perturbation expansion shown in 
Fig. 1(a) is easily demonstrated to be equal to: 

G N = (2N - 1)!! = J>™-V[ y (2) 

which is just the total number of ways in which we can connect 2N vertices by N impurity 
lines. The analogous number of graphs for in the expansion of Fig. 1(b) is much more 
difficult to find and the exact answer for this problem is, as far as we know, absent in the 
literature. Only in Ref. a simple inequality was shown to hold: 

(2N- 1)!! > Z N > (2JV-3)!! (3) 

This inequality gives only a very rough estimate of E^r. As we shall show below this problem 
can be solved exactly. This follows immediately from an exact solution for the electron in the 
random potential V(r) = V, where the value of V is independent of the spatial coordinate r, 
though Gaussian with distribution width < V 2 >= W 2 . Naturally, in this case the diagram 
technique is of the standard form like in Fig.l, and each impurity line transfers just the zero 
value of momentum, i.e. it is associated (in momentum space) with correlator of the form 
(2 / jr) d W 2 5(q L ) {d - is spatial dimensionality) |7]f|. All contributions of the same order in 
perturbation expansion of Fig. 1(a) are equal and the series for the Green's function can be 
expressed as j7|: 

G(E,p) = G (E,p) (l + f) (2N - l)\\G 2 N (E,p)W 2N ) (4) 
I jv=i J 

Then using the representation: 



oo 

,2 

2N-2, 



{2N - 1)!! = / dtt 2N - 2 e-- (5) 

v 2n ' 



the series (£|) is summed directly and we obtain: [] 



Mathematically it is equivalent to Borel summation. 
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W \WG (E,p) 

where we introduced the function: 



G(E,P) = -4* f , * J (6) 



-OO 



Consider the self-energy part corresponding to the Green's function (^|). Addition of 
a new impurity line to any diagram in this problem leads just to the appearance of an 
additional factor W 2 G 2 , and the self-energy part, defined by the expansion shown in Fig. 1(b) 
can be written as: 

£ = Q(W 2 G 2 )W 2 G (8) 

where Q(x) is some function. We shall see, that this function is the generating function for 
the number of "skeleton" graphs for the self-energy part, i.e. the coefficients of its Taylor 
expansion determine the required numbers Ejv- 

The Dyson's equation for the problem under consideration is written as: 

G = G + G Q m = G (l + Q(W 2 G 2 )W 2 G 2 ) (9) 

Introducing z = (WGo)' 1 and y = W 2 G 2 , from Eqs. (^) and (P) we obtain the following 
parametric representation for Q(y): 

l + yQ(y) = z^(z) =z^j 

y = V 2 (z) (10) 

This representation of Q is rather inconvenient. Below we show that this function obeys 
certain differential equation. It is easily seen that previously introduced function ^(z), 
satisfies the usual dispersion relation:^ 

„ T , s 1 f , Im^(t) 1 t / \ It 2 . . 

Rety(z) = - / dt — ; -ImV(t) = T^=e~~ 11 

7T J t — Z IX V27T 



-oo 



2 The sign of imaginary part is defined as we consider either retarded or advanced Green's function. 
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from which it follows immediately that *&(z) satisfies the differential equation: 

dm 



with the initial condition: 



(12) 

dz 



*(* = ±»*0) = (13) 
Differentiating the first equation in ([!(]) over y, we obtain: 

^H^T^- 1 } (14) 

Differentiating the second-equation in C [3~0|) over z and using (|i2|) we obtain: 

^ = 2*(^)^ = 2¥(z)(l - = -2yf g(y) (15) 

Comparing (|14D and (|15|) we obtain the following non-linear differential equation for Q(y): 

^Sr^-o-'fr) +*?<»)} (w) 

Using (|ig) and (UJ) we get y = ^ 2 (^ =±i0 = -f , so that 



= I ( 17 ) 

2=±r0 



2' 2/ 

which gives the initial condition for Eq. ([l6|) . Note that the condition Q(0) = 1, which 
obviously follows from diagram expansion for E, is a special one for Eq.(^) and can not be 
used as an initial condition. 

Eq.fll6|) can be expressed in more convenient form: 

Q(y)=l+y±yQ*(y) (18) 
We are interested in Taylor expansion for Q(y): 

Q(y) = E a ny n (19) 

n=0 

As the number of "skeleton" diagrams of the iV-th order for the self-energy part is just the 
coefficient of W 2N in the expansion of £ over the powers of W 2 , it is easily seen that Eq.(^) 
gives the required value of E^r as: 



Sat = a N -i (20) 

This means that the function Q(y) is the generating function for combinatorial factors Sjy. 

Substitution of Eq.(|19D into Eq.(|18D leads to the following recurrence relation for the 
coefficients a n : 

n-l 

Q"n = ^ ] ^m^n-l-m (21) 
m=0 

where ao = 1. From ao = 1 it follows that Q(0) = 1. It is due to this fact this point is 



special - equation Q(0) = 1 is satisfied for any initial conditions, for which Eq . (JIS) has a 
solution. 

From Eq.(|2lD it is easy to find the values of a n for small n, the appropriate results are 
presented in Table 1. 

The knowledge of combinatorics for the self-energy part allows to find also the combi- 
natorics for the two-particle Green's function - both for the full vertex-part T and for the 
irreducible vertex U. The appropriate diagram representations of these vertices is shown in 
Fig. 2. The self-energy part S is connected with the vertex-part T by the equation shown 
graphically in Fig.3. For the problem with zero transferred momentum 0,^] this equation 
has the following form: 

E = W 2 G(1 + G 2 r) (22) 

Thus, for the number of diagrams in the N-th order of the full vertex T^, we obtain imme- 
diately: 

Tat = Sjv +1 = ajy (23) 

In this sense Q(y) is also the generating function for the number of diagrams for the full 
vertex-part. 

The number of diagrams of the N-th order for the irreducible vertex-part Un can be 
easily obtained if we note, that the cut of any of 2N — 1 internal Green's function lines in 
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the diagram for the self-energy part of iV-th order produces the appropriate diagram for the 
iV-th order contribution to U (Fig. 4). Thus: 

U N = (2N - 1)Ejv = (2N - l)ajv-i (24) 



In Appendix A we once again derive the differential equation ( IB ) for the generating 
function Q(y), using only the Bethe-Salpeter equation, connecting U and T, as well as the 
Ward- type identity, with no use of the explicit solution for the Green's function 

III. ASYMTOTICS FOR THE NUMBER OF DIAGRAMS IN THE LIMIT OF 

LARGE N. 

In the limit of large orders of perturbation theory N ^> 1 the use of the recurrence 
relation fl2~l"|) becomes inconvenient due to the factorial growth of the number of diagrams 
||. At the same time the fact of this factorial growth itself can be used for significant 
simplification of the problem. Let us rewrite (|21"D as: 



a n = 2na a„-i + 2na 1 a n - 2 + 2na 2 a n - 3 + ■ ■ ■ (25) 

where ao = 1, a% = 1, a 2 = 4. It is natural to assume that in the limit of large n we 
have a n ~ (2n + P)a n _i, then a n _ 2 ~ 2 n-2+/3 e ^ 0. Substituting these expressions into (|25|) 
immediately leads to /3 = 1 and 

a n = {2n + l + 0{-))a n _ l (26) 
n 

This means that in the limit of large n we have a n ~ (2n + 1)!!. Let us define b n as: 

6 -(2^ (2?) 
Substituting (^) into (pl|), we obtain the recurrence relation for 6 n : 

^ (2m+l)!!(2n-2m-l)!!, , 
= n E " (2n + l)!! ~ h ^~^ ^ 



S 



with bo = 1. In the limit of large n, taking into account b\ = |, b 2 = jg, and limiting 

ourselves to the accuracy of the order of b/n 2 (where b ~ b n ~ b n -i ~ b ), we get: 

A6 n = 6 n -6 n _ 1 = ^ + 0(A) (29) 
Thus, in the limit of large n we can write down the following differential equation for b n : 

M» = lhL + o(±-) (30) 
dn An 2 n 3 

from which it follows immediately: 



bn = b ■ exp ( -\- + 0(1)) =b(l-j- + 0(1) \ (31) 
4 n / I An n l 



It is natural, that this analysis can not provide us with the value of the constant b = 
lim n _ +00 6 n . Numerical study of b n using the recurrence relation (|28|) completely supports 
the dependence defined by Eq . (|3l|) (Cf.Fig.5) and leads to the value of b — - — 0.36787944 • • • 
(calculations were done up to n = 5000, which guarantees the claimed accuracy). We do 
not know any analytical way to obtain this rather curious result. 

Finally, the asymptotic expressions for the number of diagrams of different types for 
large N have the following form: Q 

Z N = a N ^ = b N ^(2N - 1)!! = ~ {l - ~1 + 0(1^)} (2N - 1)!! = 

1_ (i_£l + O(l-) )2 JV r(iV + -) (32) 
/ne I AN V iV 2; J V 2> y 1 



r N = a N = l{l-l± + 0(±)}(2N + l)U 



' f ~ + 0(-L)\2^T(N + 3 -) (33) 



0Fe I AN y N 2 ') v 2 



17* = {2N - l)a N ^ = i {l - ^1 + 0(1^)} (2N - l)(2N - 1)!! = 



3 Asymptotic dependence of the type of Eq.(p2|) S^r c2 A, T(iV + /?) was obtained in Ref. 
However, the values of the coefficients c and (3 were not determined. 



It is interesting to note that: 




In Table 1 we present the summary of the main results for the number of diagrams of 
different types. 



IV. ELECTRON IN THE GAUSSIAN RANDOM FIELD WITH THE 
"WHITE-NOISE" CORRELATOR. 

As an example of the practical use of the results obtained above, let us consider the 
problem of an electron in the Gaussian random field with "white-noise" correlator, when 
the expression associated with impurity interaction line is |[|,|2|||: 

^(Pl, P2, P3, P4) = W 2 5(p! - p 2 + P3 - P4) (37) 

where W 2 = pV 2 , p - density of impurities, V - Born scattering amplitude of the point-like 
impurity. It is well known that the basic difficulties in this problem appear in the energy 
region, defined by the condition M: 

\E\ ~ j(E) \E\ ~ E sc (38) 

where j(E) = irpV 2 N(E) - is Born scattering rate (N(E) - density of states at the energy 
E), E sc ~ m'^d^pV 2 )^^ - characteristic size of the "critical" region around the band-edge, 
where the strong scattering appears. These difficulties are mainly due to the impossibility to 
sum any kind of dominating series of Feynman diagrams, analogous to that being summed 
in the limit of weak scattering E ^> j(E), E ^> E sc [} In fact all the diagrams for the 



4 In this case diagrams with noncrossing impurity lines dominate and we have only to take into 
account the first diagram in Fig. 1(b) 
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self-energy part become of the same order of magnitude in the energy region \E\ ~ E sc and 
should be taken into account. 

Perturbation expansion for the self-energy part in terms of "skeleton" graphs is shown 
in Fig. 1(b). All the graphs of 3-rd order in this expansion are in fact equal to each other 
(Appendix B). Despite the fact that this property is lost already in the next order, it seems 
reasonable to formulate an approximation assuming the equality of all diagrams of this type 
in every order of perturbation theory. Apparently this approximation can be good enough 
especially in the "critical region" \E\ ~ E sc where all contributions are of the same order 
of magnitude. Let us take as a "basic" graph in any given order the "maximally crossed" 
diagram like shown in Fig.6(a). For systems invariant to time reversal this graph can be 
transformed into a "ladder" -like, as shown in Fig.6(b). Then the full series for the self-energy 
part in our approximation can be expressed like: 



n=l 



E(p) = £ W 2 Z n £ £ [W 2 G( P1 + P2 + p)G(- P2 )] B_1 G(p 



Pl P2 



pi 



W 2 £G(pi-p 2 + p)G(p 2 ) 



P2 



G(pi 



(39) 



where we used notations defined in (|19"D and (|0j), as well as the property G(p) = G(— p), 
valid for the isotropic system. Accordingly we obtain closed equation for the averaged one- 
particle Green's function: 



G- 1 (p) = G 1 (p)-W 2 J2Q 



^ 2 £G(Pi-q)G( Pl ) 



pi 



G(p + q) 



(40) 



where G^(p) = E - ^. All the non-trivial part of the problem is contained now in our 
generating function Q(y), which defines the "kernel" of the complicated nonlinear integral 
equation (50). Naturally, if we limit ourselves to the first term in the expansion (TO), we 
get Q = 1 and Eq.(^UJ) reduces to the standard sum of "noncrossing" diagrams The 
obvious advantage of Eq.(ID) in comparison with standard approach (T]||, based upon the 
summation of dominating diagrams (e.g. accounting only the first graph of Fig. 1(b)) is the 
formal account of all diagrams, which is made, however, in the approximation, assuming the 
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equality of all the "skeleton" graphs for the self-energy part in a given order of perturbation 
theory 

Equation ([10) is quite complicated nonlinear integral equation and can not be solved in 
general case, more so due to the fact that we do not know the general form of Q(y) (which 
enters Eq.([H]) as a function of a complex argument). Below we shall limit ourselves with 
some qualitative analysis of Eq.(^). Let us write Eq. fljQf) in the following compact form: 



G-\p) = G \p) -W 2 Q [W 2 G®G] ®G (41) 

where we have introduced the generalized product (convolution) of functions as: 

F®$ = ^F(p-q)$(p) (42) 
p 

Let us return to Eqs.(|i~0|), defining Q parametrically. The second equation in ( JTOj ) can be 
written now as: 

G ® G = ^^ Z ) ( 43 ) 

We have seen that in the problem with zero transferred momentum z = W~ 1 Gq 1 . Consider 
now the limit of W — > in (f43|). Then the left-hand side of fl4~3| ) reduces to Go <8> Go, while 
in the right-hand side we can assume z ~ W~ x in analogy with the problem with zero 
transferred momentum and use the easily established asymptotics ^(z) ~ =■ for \z\ ^> 1. 
Here we are slightly inaccurate, since the exact form of *&(z) is: 

9(z) = R{z)=Fi^e-£ (44) 

where for R(z) we have asymptotic expansion: 



R(z) = e 2 / e 2 dt = - + — + — H (— - < argz < -) (45) 



o 



We are using the asymptotics *&(z) ~ -, which is not rigorous, however, the results obtained 
below using this approximation are confirmed by more accurate, but rather long analysis. 
Thus, in the limit of W — > Eq. (|43|) reduces to: 
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G ® G 



1 

W 2 z 2 



1 + 0(W 2 



Wy/Go ® Go 

Accordingly, in the limit of — > instead of (fi|) we can write: 



(46) 



G®G 



1 



1 



(47) 



IV 2 ^ \Wy/G ®G 0/ 

Consider the energy region E < 0, where the fluctuation "tail" in the density of states 
appears [P|,pHf|. In this case from (^) we obviously have z G Re. With the help of (|44| ) and 
(pEBj) we obtain from (pnj): 

2 ^ ^ r ii 



G ® G « G ® G - %—^\jGo® G exp j- 



If V2 V u ~ ' I 2^ 2 G ®G . 
where the second term, as we shall see now, produces the fluctuation "tail" in the density 
of states, using J2qG £§) G = D P E q G(p — q)G(p) = fe p G(p)) , we immediately obtain 
from (|48|) the density of states in the form: 

■} 



N(E) 



1 



2W 2 G ®G Q 



7T 



EpGo(^,p) 



(49) 



Now everything is defined by the concrete form of Go <S> Go for different spatial dimensions. 
The denominator of (49) is easily calculated in general 

l 



pi; 

Sd J dpp 
o 



d-1 



\E\ + f- 



(50) 



where Sd = 2~^~ 1 ^7r d / 2 r( ,J^ 2 ^ while the value of po - represents the cutoff of the order of 
inverse interatomic distance |J , which is necessary for d > 2 (for d— I'm ([5(]) we can extend 
integration up to infinity). Now E in Eq.(|50"D denotes the renormalized energy, calculated 
with respect to the band-edge, determined in "one- loop" approximation 0, which takes into 
account the infinite (in the limit of po — > oo (for d > 2)) shift of this edge. The integral in 
(pOj) is easily calculated for any value of d. 

In one-dimensional (d = 1) case all the integrals, contributing to (^), are calculated 
exactly. After rather tedious, but elementary, calculations we obtain (Appendix C): 



\E\I 



(51) 
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The exponential factor in ( |5TD differs from the well known exact result of Halperin QX1|] (Cf. 
also Ch.II in (TIJ) by the absence of the factor of |. Preexponential factor in (|51|) is also 
different from the exact value which is of the order of ~ |yi • However, the qualitative 
form of the "tail" in the density of states is reproduced more or less satisfactorily, despite 
rather common opinion, that this "tail" can not be obtained using the usual perturbation 
theory. 

Similar, but approximate, calculation of the density of states via Eq.(^) fro d = 3 
(Appendix D) gives: 

N(E) = ^{2m)-^±-exp j ) (52) 

2 

where we have introduced the cutoff energy E Q = J^- ||. Here again the exponential factor 
coincides (up to a constant) with known results of non-perturbative instanton approach 
PJT2]-pH[ , while the preexponential in ( |52"D does not agree with any of the variants obtained 



in these papers. However, in general, the result of Eq.(^) is again satisfactory enough, 
taking into account the approximate nature of our Eq. (^0|) . f] 

Especially interesting is the analysis of Eq . (f40|) in the "strong coupling" region ||, which 
is defined by (0), i.e. in the vicinity of the initial band-edge, where the transition from 
extended to localized states takes place. In this region it is reasonable to assume, that the 
approximation of equal "skeleton" graphs contributions to the self-energy in the given order 
of perturbation theory can be rather good, simply due to the known fact that they are of the 
same order of magnitude. Strong condition of the type of Eq.(|38|) is obviously equivalent to 
the limit of W — > oo. In this limit, in "zero-order" approximation we can neglect the first 
term in the r.h.s. of Eq.(fTl|) and write: 



5 For d > 4 the use of asymptotic expression of Eq.p^) and statistical analysis of Ref. || allows to 
find the correct power of W^ 1 in the preexponential factor of the density of states. In this case our 
approximation is equivalent to the hypothesis of stationarity of higher-order contributions used in 
Ref. ||, which is valid for d > 4. 
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G (p) = -W 2 Q W 2 G®G 



G 



(53) 



It is easy to convince oneself that this is equivalent to the limit of z = ±i0 in 
in fllUp. In this case (f|3|) reduces to: 

W 2 G ®G = ^(z = ±i0) = 

and from (IT7D we get: 



or y 



7i 

2 



Q 



W 2 G®G 



2 

7i 



Formal solution of Eq. fl54|) has the form: 



G = ±i 



(54) 



(55) 



(56) 



where N = ^ P 1 is the number of states in the band. Direct substitution of (|56f) and ( 55]) 
into (E^) shows that this equation is satisfied. Thus in the "first-order" approximation, in 



the limit of W — > oo we can write down the Green's function (|41|) as: 

1 



G(p) 



(57) 



G^fo)- jf^Ep G(p) 
which is surprisingly coincides with the result of the self-consistent Born approximation 
(equivalent to the first diagram in Fig. 1(b) or Fig. 3) only with an extra factor of -. 

Obviously Eq.([57]) leads to the density of states of the Born approximation Nq(E), which 
practically coincides for d = 3 with the density of states of free electrons (with the account 
of the shift of the band-edge in one- loop approximation). In Fig. 7 we compare the results 
following from Eq.fl57|) for the density of states in one- dimensional (d = 1) system with the 
exact result of Halperin (Tl| , which demonstrates rather satisfactory agreement in the region 
of "strong coupling" \E\ < E sc ~ msWs. The width of this region grows with the growth 
of W. Note that while the "tail" of the density of states is suppressed with the growth of 
W (Cf.(pT])), the crossover region, where \E\ ~ E sc , becomes wider. 

It is possible, that Eq.(|57|) provides the qualitative justification of the use of the simplest 
Born-like approximation for the Green's function in the approaches similar to the self- 
consistent theory of localization f9||l5| - the mobility edge appears in the "strong coupling" 
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region \E\ ~ E sc (|38|), where the approximations leading to (|57|) becomes rather satisfactory 
and the Green's function really takes the simple Born-like form. 

These results demonstrate the effectiveness of the use of diagram combinatorics in for- 
mulation of new approximations for realistic physical problems. 

This work is partially supported by the Russian Foundation for Basic Research un- 
der the grant No. 96-02-16065, as well as from the Project No.IX.l of the State Program 
"Statistical Physics" of the Russian Ministry of Science. The authors are grateful to Dr. 
A.I.Posazhennikova for the help with numerical calculations. 
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APPENDIX A: 



Let us derive Eq.(jT8|) for generating function Q{y) without the use of the explicit form 
of one-particle Green's function @. In the problem with zero transferred momentum the 
Bethe-Salpeter equation shown in Fug. 2(c) reduces to: 



so that 



r = u + ug 2 t 



u 



(Al) 



(A2) 



1 - UG 2 

Using ( |A2| ) and (^) we obtain an equation, connecting the self-energy part with irreducible 
vertex U : 

W 2 G 



Use now the "Ward identity" : 







dW 



1 - UG 2 



— = UG 
W 



(A3) 



(A4) 



G 



which can be derived via (§) and (21), and Eq. ( |A3| ) to obtain: 



W 1 



d 



dW 



G 



E = W 2 G + W 2 GY, 



d 



dW 



G 



E 

w 



(A5) 



Now from (^), we obtain our differential equation for Q: 

d 



Q(W 2 G 2 ) = 1 + W 2 GQ(W 2 G 2 ) 



dW 



WGQ(W 2 G 2 ) 



G 



1 + W 2 G 



d 



-W 2 G 2 Q 2 (W 2 G 2 " 



d(W 2 G 2 



which can be rewritten as: 



Q(y)=l+y±yQ 2 (y) 

dy 



(A6) 



Note, however, that from this analysis it is impossible to find the correct initial condition 
(|17|), which is closely related to (0), reflecting the causality principle. 
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APPENDIX B: 



Let us show that all diagrams of the 3-rd order, shown in Fig. 1(b) are equal. Consider 
the "maximally crossed" graph of Fig.l(b.l). Its analytical form is: 

Fig.l(b.l) = W 6 £ G( Pl )G( P2 )G( P3 )G( P ,)G(p 5 ) x 

Pl,P2,P3:P4,P5 

X *(P - Pl + P3 - P4)£(Pl - P2 + P4 - P5)<5(P2 - P3 + P5 - P) (Bl) 

In isotropic system the Green's function depends only on the absolute magnitude of the 
momentum and G(jp) = G(—p). In ( |ET| ) we can make the following change of variables of 
integration (which obviously does not change the value of integrals): pi ^ — P3; P2 — P2- 
Correspondingly: 

<*(P - Pl + P3 - P4)5(Pl - P2 + P4 - P5)5(p 2 - P3 + P5 ~ P) ~> 
-> 8(p - Pl +P3 - P4)5(P2 - P3 + P4 - Ps)£(Pl - P2 + P5 - P) (B2) 

and we obtain the contribution of diagram shown in Fig.l(b.2). Analogously, the change of 
variables p 2 — P4; P3 — P3 produces the contribution of diagram of Fig.l(b.4), while 
P3 ^ ~P5] P4 ^ — P4 gives the contribution of Fig.l(b.3). Thus, all the contributions of 
"skeleton" diagrams of the 3-rd order are equal to each other. 

Analogous change of integration variables allows to show that under any impurity in 
the given "skeleton" diagram we can actually "rotate" the electronic line with all incoming 
impurity lines, so that the new diagram will give the same contribution as an initial one. 
This symmetry produces rather wide classes of equal diagrams in any order of perturbation 
theory. 
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APPENDIX C: 



Here we present some details of calculations, leading from (f49|) to ([51]) for one-dimensional 
(d = 1) case. For E < the convolution of Green's functions, entering (^9|) is calculated by 
usual contour integration: 



G ® G 



2tt 



dp- 



1 



1 



E\ + p 2 -2pq + q 2 \E\ + p 2 J\E~\(q 2 + 4\E\) 



(CI) 



(Here and below for brevity we put 2m = 1.) Then: 



\fGo® G exp | - 

n v. 



2W 2 G (8) G c 



7T 



1 



\E\ 



2tt\E\4 



exp 



W 2 



\E\i 
W 2 



E\ 



W 



dq 



q 2 + A\E\ 



-exp 



2W 2 



+ 4|£|) 



2\E\ 



2V2n\E 



■exp 



After that by elementary calculations we get: 



EGo(p) 



dp-. 



1 



W 2 



1 



(C2) 



\E\+p 2 2 x f\E\ 

u V 

Putting all these expressions into ([19|) and restoring 2m, we obtain ([51 



(C3) 
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APPENDIX D: 



Here we present the details of calculations leading from (pE9|) to (|5^) for d = 3. Calcula- 
tions in this case are rather tedious and it is convenient to use the coordinate representation 
for the Green's functions. In the region of E < we finally obtain: 



q < 2J\E 



j_ 



(Dl) 



q > 2yJ\E\ 

In inequalities in ( pi| ) it is convenient to replace 2 by tt, which guarantees the smooth 
crossover between two limiting expressions of Eq. ( |D 1| ) in the region of intermediate q. Using 
we can write: 



\JGq® G exp | - 

n ^ 



1 



2tt 2 



dqq 2 



1 



2iy 2 G ® G . 

^\I\E\ 



exp 



TT 



\E\*{l + 0(W 2 )}exp 



W 2 
4n x f\E~\ 



OO i 



W\E\ 



In fact the main contribution here is from the first integral. Then: 



(D2) 



1 po 2 
If, p 
dp 



71^ 







E\+p 2 2vr 



1 / 7T 

Po - -\/\E\ 



Po 
2tt 2 



(D3) 



Putting all these expressions into Eq. fl49|) and restoring 2m, we obtain fl52|) . 
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Figure Captions. 



Fig.l. Diagrammatic expansion for the averaged one-particle Green's function (a) and 
self-energy part (b). Dashed line corresponds to the quadratic correlator of the random 
field. Go ~~ free particle Green's function. 

Fig. 2. Diagrammatic expansion for the full vertex-part F (a), for the irreducible 
vertex-part U (b) and Bethe-Salpeter equation, connecting T and U (c). 

Fig. 3. Equation connecting the self-energy part with the full vertex-part. 

Fig. 4. The cut of any of 2N — 1 internal Green's function lines in the "skeleton" diagram 
of iV-th order for the self-energy part produces the appropriate diagram for U. 

Fig. 5. The behavior of b n with the growth of n. Dots represent the values of b n , obtained 



from the recurrence relation (28[), the curve represents the asymptotic dependence 



e {} ~~ In)' hashed line — the asymptote ~. 



Fig. 6. (a) - "basic" diagram, used to construct the approximation for the self-energy part, 
(b) - "rotated maximally crossed" diagrams produce the "ladder" in case of time reversal 
invariant system. 

Fig. 7. Density of states in one- dimensional system for the different values of the averaged 
square of the random field — ■ \ 0.25; 2. 2; 3. 16. Continuous curves — the exact 
solution, dashed lines — self-consistent Born approximation (0). Energy is in units of E , 



density of states - in units of ^7g=> E - some arbitrary unit of energy. 



23 



+ 



+ 



a- 1 
x5 



f 

-J 1 

i 

\ - ■ 



'1 /■ 



I \ I' 



+ 



\ 1 



■l 



;. I 

I— "1 
+ 

I--! 



to 



V 



T 



4- 



i I 




2 



i i-n 



Si 




n 



0,35 
0,30 
0,25 
0,20 
0,15 
0,10 
0,05 
0,00 



